Two-parameter scaling law of the Anderson transition 
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It is shown that the Anderson transition (AT) in 3d obeys a two-parameter scaling law, derived 
from a pair of anisotropic scaling transformations, and corresponding critical exponents and scal- 
ing function calculated, using a high-precision numerical finite-size scaling study of the smallest 
Lyapunov exponent of quasi- Id systems of rectangular cross-section of L x SL atoms in the limit 
of infinite L and SL < L, for x — SL/L ranging from 1/30 to 1/4. The second parameter is x, 
and there are two singularities: apart from the two-parameter scaling describing AT for x > 0, 
corrections to scaling due to the irrelevant scaling field diverge when x — > 0, and the corresponding 
crossover length scale is also estimated. Furthermore, results suggest that the signatures of the AT 
in 3d should be present also in 2d strongly localized regime. 
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One of the long-standing open fundamental prob- 
lems of the physics of quantum-mechanical disordered 
systems is the quantitative description of the metal- 
insulator transition induced by the Anderson localiza- 
tion of electronic eigenstates [if, known as the Anderson 
or localization-delocalization transition (AT). The subse- 
quently developed scaling theory of localization of Abra- 
hams et al. |2j (STL) proposed that: (i) in 3d the transi- 
tion is a continuous phase transition with only one rele- 
vant scaling variable (which became known as the single- 
or one-parameter scaling hypothesis), (ii) the lower crit- 
ical dimension is 2 where all states are localized for arbi- 
trarily small finite disorder strength, and (iii) the transi- 
tion can be described in terms of the scaling law of the 
disorder-averaged dimensionless conductance g that de- 
pends only on g, i.e dhxg/dhxL = Pd{g)- The theory 
is based on the renormalization-group (RG) considera- 
tions of the Thouless expression for g, and an additional 
calculation showing that Pd=2{g) < 0, which was corrob- 
orated by a more detailed self-consistent study Q based 
on a resummation of the perturbation theory for weak 
disorder 0]. 

The discovery of the universal conductance fluctua- 
tions [H, [|| showed that g is not a self-averaging quantity, 
and therefore scaling of its whole distribution has to be 
studied. The STL nevertheless survived in the sense that 
one can still find a single parameter g that characterizes 
scaling properties of the whole distribution, and there is 
a compelling evidence that the scaling properties of the 
distribution itself in the critical region of the transition 
are still governed by a single-parameter Q . 

Numerical studies of the transition using the transfer- 
matrix method (TMM) L 8| suggested a possibility that 
Pd=2(g) — for a finite amount of disorder, which would 
correspond to the existence of a line of critical points for 
disorder weaker than a certain finite disorder strength, 
but subsequent studies showed that the dependence of 
localization length on the disorder strength in 2d is in a 
quantitative agreement with analytical results 

High precision numerical calculations of the critical ex- 
ponent v describing the divergence of the localization 



leng th £ at the critical point in 3d, however, give v ks 1.54 
lfj, [ll|, in sharp contrast with v = 1 obtained from the 
self-consistent analytical calculations 0, 0| • 

Theoretical breakthrough was made by Efetov 
who introduced supersymmetry to calculate disorder- 
averaged products of Green's functions, and provided a 
theoretical framework that, among other results, allows 
calculation of v beyond the self-consistent approach, al- 
though technical difficulties with the e-expansion do not 
yet permit accurate estimate of v despite the consider- 
able theoretical progress [T3|, and currently v is most 
accurately determined using TMM. 

This Letter present three main results: (1) there is 
an additional scaling parameter x in d = 3 describing 
the thickness to width ratio of long quasi-ld wires, and 
the corresponding two-parameter scaling law and criti- 
cal exponents are estimated numerically; (2) results are 
in agreement with (ii) and additionally show that there 
should be possible to see signatures of 3d transition in 
2d; (3) the description of the transition in terms of the 
/3-funetion depending only on L is incomplete due to the 
geometric nature of both x and L. 

The starting point is the Anderson model [l|: 
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where represents the impurity energy at site i. £j is 
randomly, independently and uniformly distributed in 
[— W/2, W/2); t is the hopping integral of electron (set 
to 1), (ij) denotes that the hopping takes place only be- 
tween the nearest-neighbors of the simple cubic lattice, 
and the Fermi energy is set to (band center). 

The geometry studied is that of the quasi- lei slabs of 
L X SL X M atoms with M 3> L, ratio x — SL/L, open 
boundary conditions (b.c) in SL direction and periodic 
b.c in L direction. Similar geometry of cubic samples of 
Lo x Lo x (qLo) atoms was studied in Ref. [3], where 
authors found that the critical disorder strength W c is 
approx. independent of the shape and that g becomes 
strongly suppressed for g< 1 and g> 1 which is consis- 
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tent with, respectively, the quasi-2<i and quasi-ld char- 
acter of samples in these cases. 

The scaling properties of AT in 3d are studied by the 
standard calculation of the smallest Lyapunov exponent 
7 of transfer matrices of long quasi-ld samples [1(1 flit - 
The inverse of 7 is the largest length scale in the problem, 
which is identified with the correlation length £ = 1/7. 
The usually studied quantity is the rescaled correlation 
length A defined as A(L, W) = (Lj(L, W))" 1 

The finite-size scaling analysis of A gives the scaling 
properties of 3d systems in the thermodynamic L — > 00 
limit, by considering how A(L, W) changes under the RG 
transformation R:Jjh b'L,SL b'L [HI. The corre- 
sponding scaling law, including the corrections to scaling 
due to one irrelevant field was considered in the context 
of AT first by Slevin and Othsuki 
numerically show that 



ll|, who were able to 



A(L,w) = F(L»i>{w),L v (t>{w)), w 



W-W c 



(2) 



where v > and y < are critical exponents associ- 
ated with, respectively, the relevant and irrelevant scal- 
ing fields tp and <j), and all functions are analytic. This 
is achieved by fitting numerically obtained values of A 
with the truncated expansions of F, 4>, and %p, while the 
error-bars are estimated via calculation of 95% confidence 
intervals that can be done either by the bootstrapping 
method [llj] or by a direct calculation of projections of 
the confidence region [lj|. 

In the case when there is an additional parameter x, 
we can repeated the above procedure for several values 
of x, 



A(L, 10, x) = F(L"ip(w, x), L v 4>(w, x),x), 



(3) 



where v,y,W c may in general also depend on x. The 
principal difference between Eq. @ and Eq. is that 
F, ip and <fi do not have to be necessarily analytic in x. 
Expansion in the second argument of F gives: 

A(L,w,x) = F±(L"ij)(w, x), x), A = A — AF, (4) 

oc 

AF = J2 F n( Lh ^,x),x)^(w,x) n L n y, (5) 
n=l 

where F± represents the universal part describing AT, 
and AF are corrections to scaling due to 4>. These van- 
ish for large L because y < but are important for a 
quantitative description of AT, including precise deter- 
mination of all of the relevant parameters [lj [l6| . 

Figure 1 shows the typical behavior of A(L, W) and 
the corresponding fit for small constant x and several 5L 
starting with 1. As L = SL/x increases, A(L, W) at first 
decreases uniformly in W (since for small SL system is 
close to being 2d) but with further increase of I a char- 
acteristic behavior for the continuous transition begins 
to develop, since for large SL and constant x system is 
3d. 
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FIG. 1: Values of A(L, W) for x = 1/10. The corresponding 
fits are represented with lines for 5L = 1, 6 (going from the 
slowest changing curve to the fastest). 



Table U summarizes the parameters and the results of 
the numerical simulation using the methods described 
above for several values of x. Results presented in the 
Table suggest that v is approx. independent of x and 
in agreement with Ref. [ToL 11 . Il6j |. which supports its 



universality. Similar can be said for W c , in agreement 
with Ref. 1^|, while y becomes approx. constant for x is 
1/4. 

To better understand the scaling properties of AT w.r.t 
x, I consider an additional scaling transformation, 1Z X : 
x 1 ► bx under constant L (this is equivalent to scaling 
only the thickness SL), and introduce a scaling field ip 
that depends on x but not on w: 



A = F±(L^(w,x),tp(x)), 



(6) 



where W c ,v and y are assumed to be constant. The 
scaling law of F± w.r.t x can be derived assuming that 
functions scale under 1Z X in the general way [TBI ] : 



K X :F±^ b yi F ±1 i-> b V2 tP, V 1 — ► b V3 ip : 



(7) 



Iterating TZ X a finite number of times in the standard 
manner gives the two-parameter scaling law: 



A = i a f ± (i5^(w)), 



(8) 



where a = -yi/y 3 , fi = -2/2/2/3- 

At first it seems that Eq. (J8]) cannot be correct in the 
case of AT for small x because numerical results give that 
the transition takes place at approx. const. W c , and if 
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SL 


W 


N d 


N p 


x 2 


Q 


W c 


V 


y 


1 


[4, 14] 


[15, 18] 


427 


10 


439.1 


0.2 


16.46(39, 55) 


1.60(55, 65) 


-1.39(1.16, 1.62) 


3/4 


[3, 15] 


[15, 18] 


304 


10 


282.4 


0.7 


16.52(46, 60) 


1.62(57, 66) 


-1.13(0.98, 1.29) 


1/3 


[1,4] 


[15, 18] 


427 


10 


418.2 


0.5 


16.52(42, 63) 


1.54(48,61) 


-1.00(0.90, 1.10) 


1/4 


[1,7] 


[15, 18] 


427 


10 


436.9 


0.2 


16.48(42, 54) 


1.55(51,60) 


-1.00(0.97, 1.04) 


1/7 


[1,5] 


[14.25, 18] 


380 


11 


370.0 


0.5 


16.61(47, 77) 


1.52(39, 67) 


-0.89(0.84,0.95) 


1/10 


[1,6] 


[14.5, 18] 


370 


11 


365.7 


0.4 


16.62(50, 75) 


1.52(41,66) 


-0.88(0.84,0.93) 


1/15 


[1,4] 


[14.5, 18.25] 


304 


11 


298.2 


0.4 


16.58(23, 95) 


1.62(41,87) 


-0.90(0.80,0.99) 


1/20 


[1,4] 


[14.5, 18.25] 


304 


11 


287.3 


0.6 


16.59(27, 99) 


1.57(37, 89) 


-0.90(0.80, 1.00) 


1/30 


[1,4] 


[14.5, 18.5] 


324 


11 


350.3 


0.1 


16.50(27, 80) 


1.59(51,69) 


-0.93(0.85, 1.01) 



TABLE I: Input parameters of the simulation: x, thickness SL, the number of points Nd, the number of parameters N p ; Output 
values: the minimal least-square deviation x 2 obtained, the quality of fit parameter Q; and the values of the fitting parameters: 
the critical disorder strength W c , the critical exponent v and the irrelevant exponent y. All A's are calculated with 0.1% relative 
accuracy except for the largest systems studied (x — 1/30, SL = 4), which is calculated with accuracy of 0.12%. Error-bars are 
95% confidence intervals. 



W c would remain constant for smaller x as well, Eq. ([5]) 
would imply that transition persists when x — > (re- 
gardless of the value of a), where one must find strongly 
localized states instead. 

This is resolved if we notice that AF is also divergent 
as a; — > 0, and that there is a characteristic crossover 
length-scale L c (x) such that for L > L c (x) corrections to 
scaling become small, but L c (x) — ► oo when x — > 0. The 
source of the divergence for W ps W c is in <f>(x, w = 0) = 
</>o(x), and the left panel of Fig. [2] shows that <j){x) oc x 5 
for x < 1/7, with S = -0.17(05, 29). L c (x) is determined 
from the condition L c (x) v 4>{x) ~ 1, giving 

L (x)~x~i. (9) 

A more detailed discussion of the exponent S will be car- 
ried out elsewhere since our main interest here is in the 
two-parameter scaling law Eq. © of AT. 

a is calculated from the dependence of A c on a;, where 
A c = A(W = W c ). The right panel of Fig. H shows that 
a = 0.89(84,95). Numerical verification of Eq. © is 
done by a rescaling of the argument of A/A c for each x 
individually by a factor 6{x) chosen such that all A/A c 
collapse, and Fig.[3]shows the result. Figure Q] shows that 
0(x) oc a/ 11 , /ii = 2.75(2.28,3.16) and ip(x) oc a^ 2 ,^ = 
1.34(22,46), which gives fjt = m + /z 2 = 4.1 ± 0.6. 

The largest similarity with the one-parameter scaling 
Eq. {2]| is achieved for a = 1, fx = when Eq. {SJ, 
expressed in terms of £, becomes 

Z{L,w,x) ~ (xL)T ± ((xL)^j(w)). (10) 

x influences the divergence of the correlation length of 
the infinite system £,co(w,x) ~ (\w\/x)~v for L > L c {x). 
The numerical results are not incompatible with a = 1, 
but nonetheless strongly suggest /i ^ a/ v. 

The two-parameter scaling was obtained for samples of 
small SL, including the single-layered case (Table[T|. This 
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FIG. 2: Dependence of 4>q{x) (left panel) and A c (a;) (right 
panel) on x. The solid line in both panels is the least-square 
linear fit. Error-bars are 95% confidence intervals. 



allows one to follow outflow RG trajectories from L — > oo 
back to the single-layered finite L case, and therefore pro- 
vides an analytical mapping of the phase diagram of 3d 
systems onto 2d systems, suggesting a possibility of sig- 
natures of 3d AT also in 2d strongly-localized regime. 
What exactly those signature are, including a possibil- 
ity of their experimental observation in the mesoscopic 
regime, will be discussed elsewhere (l7| . 

Although the study carried here shows the possibility 
of AT in arbitrarily thin 3d films (which should be of 
size L > L c (x) in order to reduce the large corrections to 
scaling due to the dimensional crossover) in the strong- 
disorder regime, the critical conductance is strongly sup- 
pressed due to the x a prefactor. 

The RG approach proposed seems to be a rather accu- 
rate description of AT and allows the nontrivial extension 
of the one-parameter STL to two parameters. It remains 
to be seen whether such an approach can be applicable 
to other problems in physics, for instance in addition to 
the dimensional regularization. 
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FIG. 3: Collapse of A(L, W, x) onto the universal function 
J-±, for 2109 points from Table U corresponding to x < 1/4. 
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FIG. 4: Dependence of the rescaling factor 8(x) (left panel) 
and the relevant scaling field ip (right panel) on x. The solid 
fine in both panels is the least-square linear fit. Error-bars in 
the right panel are 95% confidence intervals. 
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